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Intermediate evolution using SNIa, and BAO 

Victor H. Cardenas and O. Herrera 


Abstract We study the intermediate evolution model 
and show that, compared with the recent study of a 
power-law evolution, the intermediate evolution is a 
better description of the low-redshift regime supported 
by observations from type la supernovae and BAO. We 
found also that recent data suggest that the intermedi¬ 
ate evolution is as good a fit to this redshift range as 
the ACDM model. 
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1 Introduction 


The current cosmological paradigm was set primar¬ 
ily after discovering that type la supernovae (SNIa), 
a well-studied standard candle, are dimmer than ex¬ 
pected in the context of t he standard cosmologica l 
model ( Riess et al. {1998}; Perlmutter et al. ('19991 1. 
The simplest way to describe this observation in the 
context of the standard big-bang picture is by intro¬ 
ducing a cosmological constant term. A, into Einstein’s 
field equations. 

In this setup the cosmological constant A drives the 
current accelerated expansion of the universe, lead¬ 
ing to the successful ACDM model, the simplest one 
that fits a varied set of observational data such as 
SNIa, measurements of the baryon acoustic oscillations 
(BAO), information from the cosmic microwave back- 
CTOund radiation fCMBRl, growth of structur e, etc 
( Weinberg et al. ('20131 : Perivolaropoulos 120101 1. 


However, this model looks unnatural in at least two 
ways; (i) there is no clue about the physical mechanism 
to get the current value for Da, and (ii) this model 
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tells US we live in a very special epoch where the cos¬ 
mological constant contribution Da is of the same or¬ 
der of magnitude as the dark matter contribution Dm, 
a highly improbable fact considering the dark matter 
contribution decreases as a“^, with a{t) the scale fac¬ 
tor. Meanwhile the cosmological constant contribution 
is and always has been a constant. 

A plausible alternative is to explore deviations from 
the ACDM model, assuming a variable A. Dark en¬ 
ergy (DE) is the name of the unknown component re¬ 
sponsible for the current acce lerated expansion of the 
universe (jErieman et al.ll2008h . In its simplest form. 


this can be described by a fluid with constant equa¬ 
tion of state (EoS) parameter w = — 1 corresponding to 
a cosmological constant. There are also models where 
a scalar field drives t he cosmic ac celeration, so-called 
quintessence models (IMartinl 120081) . as well as models 
where through a modification of the gravity secto r the 
cosmic acceleration is described ( Tsniikawa 2010l) . 

This honest approach focuses only on problem (ii), 
mentioned in the previous paragraph, but it is almost 
completely disconnected from problem (i). Most of the 
work in cosmology since the discovery of the acceler¬ 
ated expansion of the universe has been of this type. 
The efforts are focused on characterizing this new com¬ 
ponent, dark energy (DE), somehow disconnected from 
first principles. 

Many scientists think that a real understanding of 
the DE problem, focusing this time on problem (i), i.e., 
the o ld well-known c o smological constant (C.C .) prob¬ 
lem ( Weinberg 1989t Peebles and Ratra 20031) . would 
certainly reveal a door to new physics. 

Among the many attempts to tackle the cosmological 
constant problem, one of the most interesting is where 
a suitable relaxation mech anism dr i ves th e C.C. to a 
small value. Eor example, iDolgovi (|l997[) describes a 
model where a scalar field is coupled to the curvature of 
the space-time, such that their contribution to the en- 
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ergy density cancels the vacuu m energy. Similar to this 


work are those proposed by iHebecker and Wetterich 


( 2000l l. where an additional scalar field is introduced 


the G auss-Bonnet (GB) term (IBoulware fc Deseiill985 . 
1986ll . GB terms arise naturally as the leading or¬ 


which to date remains a C.G. wh ich asymptotically van¬ 
ishes (see Padmanabhajj ( 20031) for a revi ew). 


der of t he expansion to the low- energy string effective 
action ( Koivisto fc Mota 2007allb l. This kind of the¬ 


ory has been applied to a possible resolution of the 


In a recent paper (|Dolgov et al.l 1201411 the authors initial singularity problem (lAntonia dis et al Jll994r). to 


modified gravity theories at low redshift - and directly 

1993; Kanti et al. 1996: Ghen et al. 2007). and acceler- 

studied the constraints the supernova data imposed on 

ated cosmological solutions (Noiiri & Odintsov 2005: 

the parameter B. They found, surprisingly, using both 

Noiiri et al. 20051: Coenola et al.l 2006: Barrow et al. 

data from the Union 2.1 (Suzuki et al. 

2012|l and the 

20061 del CamDol|2014). among others. 


recently released jo int light-curve analysis (JLA) set 


(|Betoule et al.ll2014l) . that SNIa data suggest a /3 ~ 3/2 


assuming a fiat universe. Although successful in pro¬ 
ducing a good fit to the SNIa data, the model is un¬ 
able to describe a transition from a decelerated to an 
accelerated phase. In fact, by using a scale factor 
a{t) ~ this implies a constant deceleration parame¬ 
ter g = (1 — /3)//3. This transition is a key ingredient 
that any model must satisfy if it claims to describe the 
recent (low redshift z < 1) evolution of the universe, 
which is essentially the same redshift range as the su¬ 
pernova data span. 

It is the purpose of this paper to explore a slightly 
different scale factor evolution - the so-called intermedi¬ 
ate evolution - that interpolates between a power-law 
evolution and an exponential one, that is able to de¬ 
scribe the transition from the decelerated phase to an 
accelerated one. We use th e latest data yts o f SNIa 
(some of which was used in ( Dolgov et al. 2014h l. and 
compare our results. 

In the next section we introduce the intermediate 
evolution. In section III we describe the results of our 
study and the comparison with previous works. We end 
with a discussion section. 


2 The intermediate evolution 

The intermediate evolution was introduced first in the 
context of inflationary models as an exact solution for 
a scalar field p otential of the type V((t>) cc <t)~ 


. p otential or tne type v (0) oc 0 y ^ 
(|Barrowlll99u Barrow fc Saich 1990l : Muslinov 1990ll . 
where / is a free parameter with range 0 < / < 1. A 
potential of this form, in the context of the slow-roll 
approximation, gives a Harrison-Zel’dovich spectrum 
of density perturbations with an exact scale-invariant 
spectral index i.e., Ug = 1. 

The motivation to study an intermediate inflation¬ 
ary model comes from string/M theory. This the¬ 
ory suggests that in order to have a ghost-free ac¬ 
tion, high order curvature invariant corrections to 
the Einstein-Hilbert action must be proportional to 


Particularly interesting to this work is the finding 
that, using a GB interaction with a scalar field (j), it is 
possible to describe a DE model leading to a solution 
of th e form a{t) = oq exp [(2/(Kn)) dSanval 2007 . 
200911 . Here, K = 8ttG and n is an arbitrary constant. 
Actually, this is exactly a particular case of the inter¬ 
mediate evolution, in which the scale factor evolves as 


a(t) = bexp 


( 1 ) 


where A is a positive constant and / was introduced 
above. Thus, the expansion of the universe is slower 
than standard de Sitter inflation (a(t) = exp (Ht)), but 
faster than power law inflation (a(t) = tP;p > 1). 

Using © the Hubble function is 




and the deceleration parameter is 
1 -/ 1 


aa 

9 = -TT = -1 


/ At/' 


( 2 ) 


( 3 ) 


This last result shows immediately that this kind of 
evolution makes it possible to describe the transition 
from a decelerated to an accelerated expansion phase. 
In fact, for small t the second term dominates, which for 
/ < 1 gives a decelerating universe evolution, and for 
late (larger) times, the first factor starts to dominate, 
leading to an accelerated evolution. 

Erom a purely theoretical point of view, most of 
the efforts made using the GB scenario to tackle 
the DE problem, have focused on obtaining solutions 
which look similar to the AGDM model at the back¬ 


ground level (see (jGranda fc Jimenea 120141 1 and refer- 
ences therein). A ver y interesting result was found in 
( De Felice et al. |[2^, where by studying the so-called 
general GB gravity, the authors show in the linear per¬ 
turbations regime that there is an instability during the 
radiation and matter domination epoch. This modified 
GB gravity is equivalent to the term GB coupled to a 
scalar field, but without a kinetic term. There remains 
the uncertainty as to whether this instability appears 
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once we turn on the kinetic term, as is the case of the 
model from which the intermediate evolution emerged. 
Although this instability does not spoil out the power 
of the model to describe DE, as the authors mentioned 
in the text, it is clear that this issue must be delved into 
more de eply and also in a more general setup (see for 
example de la Cruz-Dombriz fc Saez-Gomez ( 2012ll b 

More interesting for the present work are the re- 
sults from a direct test a gainst observational data. In 
( Koivisto fc Mota 2007all the authors use data from 
solar system, type la supernova, cosmic background 
radiation, large-scale structure and nucleosynthesis. 
They found some tension with nucleosynthesis and the 
baryon acoustic scale. However, their results are based 
on working in the special case of an exponential poten¬ 
tial for the scalar field. The intermediate evolu¬ 

tion exists for a more elaborate scalar field potential, 
and this fact makes a study of this type of evolution 
imperative. 

To test the model against the observation, it is useful 
to write down all the formulae in terms of the redshift, 
z = —1 + a{to)/a{t). We also impose that a(to) = 1, 
given the relation 


In b{l + z) = —At^. 


(4) 


From this equation we can rewrite the deceleration pa¬ 
rameter m as 


qiz) = -1 + 


/-I 


1 


/ ln6(l + z)’ 
and the Hubble function can be written as 


Hiz) = Ho 


1 -b 


ln(l -I- , 

hTT 


(5) 


( 6 ) 


Usually the ratio E{z) = H{z)/Ho is what we need to 
perform the statistical analysis. It is clear that the free 
parameters of the model are: h, f and b. The original 
parameter A in dm is related to b through the age of 
the universe to, by In6 = —At^. 

Because this model shows a transition from a decel¬ 
erated phase to an accelerated one, we choose to replace 
the b parameter by the cross redshift Zc, at which q = 0. 
From ([U we get 


ln6 = - ln(l + Zc). 


(7) 


In the next section we shall use (0, written in 
terms of the free parameters h, f and Zc, to test the 
model against SNIa data, using both the Lick Ob- 
servatory Super r iova Search (LOSS) compilation set 
( GaneshalingamI 1201311 comprised of 586 SNIa, and 


the J oint Ligh-curve Analysis (JLA) set (IBetoule et al 
2014 1 with 740 points, the largest compilation so far. 

We also a dd baryon acoustic o scillation (BAO) data 
points from (IBeutler etahl 1201 il l . using the approach 
considered by lxia et*^ ( 2012l l and Dolgoy et al.l ( 2014l l .| 
From the quoted references it is found that: 

d{z) = (0.335±0.016,0.576±0.022,1.539±0.039), (8) 

for z = (0.106,0.2,0.57) respectiyely, where d{z) = 
Dv{z)/Dv{z = 0.35) with 


Dv{z) = 


{l + zfDliz) 


cz 

m7) 


1/3 


(9) 


where Da{z) is the angular diameter distance defined 

by 


Da{z) = 


dz' 


Hoil + z)Jo F(z')’ 
yalid for flat space. 

3 Testing the model 


( 10 ) 


In th is section we use the Loss compilation set ( Ganeshalingam^ 
2 OI 3 I 1 and the JLA set ( Betoule et al. 2014 1. together 
with the BAO points to constrain the free parameters 
for the intermediate eyolution model. 


In the case of the Loss compiled set (jGaneshalingam 


2013 ), the fitting is done by minimizing the function, 
{pi - bi‘th{zi)Y 


= E' 


( 11 ) 


where Hi, Oi are the obseryational yalues of the distance 
modulus and its errors, while ^J^thizi) is the theoretical 
yalue of the distance modulus eyaluated at the obseryed 
redshift Zi, 


Hth{z) = m - M = 5 logiQ DLiz) + 25, 
where 

is the luminosity distance in f lat space. 


( 12 ) 


(13) 


In the case of the JLA set (|Betoule et al.ll2014ll . the 
function we minimize is 


X — (m f^th) C {fj, 


(14) 


wh ere C corresponds to the coyariance matrix released 
in (IBetoule et al.ll2014ll . and the distance modulus is 
assumed to take the form 


= m — M + aX — 7 Y, 


(15) 
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where m is the maximum apparent magnitude in the 
rest frame of the B band, X is related to the time 
stretching of the light-curves, and Y corrects the color 
at maximum brightness. In the general case, cosmology 
is r estricted together w ith the parameters M, X and Y. 
In ( Betoule et al.l[2014 ) is also given a compressed form 
of the data, where only M is left as a free parameter. 


3.1 ACDM model 


In order to have a reference model to compare with, 
we start studying the performanc e of the ACDM mode l 
against the JLA SNIa data set (jBetoule et al.l 120141 1. 
We assume a ACDM model with arbitrary curvature to 
ensure three free parameters in the fitting process; h, 
Dm and Da. The Hubble parameter for this model is 



f 


Fig. 2 We plot the contour at la and 2a for the parame¬ 
ters / and Zc after marginalization o f h. Here we have used 
data from the Loss compilation set (lGaneshalingamll2013l l 
and the BAO points discussed in the text. 


E (z) - Dm(l + z) z Daz(2 + z) + il + z) . (16) analysis they performed. After marginal 


The best fit gives Xmin = 32.76. After marginalizing 
in h, the following best fit parameters are obtained: 
Dm = 0.19 ± 0.11, Da = 0.56 ± 0.17, whic h agree with 
the values reported in ( Betoule et al.ll2ni4l l and in par¬ 
ticular with Figure 15 in that paper. 





Fig. 1 We plot the contour at la and 2a for the param¬ 
eters Dm and Da after marginal ization of h. Here w e have 
used the compressed data from (iBetoule et al.|[2^014l L 


3.2 Intermediate model 

In order to compar e with the results of the paper by 
Dolgov et al.l (|2014ll . here we describe the results using 
both SNIa data sets, together the BAO points. 

The best ht values of the parameters - using the Loss 
set together with the BAO points - are: / = 0.48±0.06, 
Zc = 0.63 ±0.07 and h = 0.675 ±0.004. Actually, these 
values remains essentially the same with and without 
taking the three BAO points into consid eration. This 
fact is also mentioned by the authors of (IDolgov et al 


ization of the parameter h, we plot the confidence limits 
of the parame ters / and Zr. i n Figl2] 


We follow iDolgov et al.l (|2014l l and use the com¬ 


pressed form of the JLA likelihood where only M is 
taken as a free parameter that is constrained together 
with the cosmological parameters. This procedure is 
equivalent to what we already did in the previous case 
using the Loss compiled data set. In that case, we con¬ 
sider h as a free parameter, and after the analysis was 
performed, it was marginalized. This is the well-know 
degeneracy between the absolute magnitude and the 
Hubble constant. 


Using the JLA data set (jBetoule et al.N2014ll alone, 
the best fit parameters are: / = 0.50 ± 0.28, Zc = 

1.2 ± 1.3 and M = 43.168 ± 0.024. Notice the large 
uncertainty in the determination of the cross redshift 
Zc- After marginalization of the nuisance parameter 
M, the confidence contour of the parameters / and Zc 
are displayed in Fig. ([3]). 

After adding the BAO points, the result gives small 
changes compared with the previous case with SNIa 
alone, with the following best-fit parameters: / = 
0.52 ± 0.21, Zc = 0.93 ± 0.71 and M = 43.175 ± 0.024. 
After marginalization of the nuisance parameter M we 
get the conhdence contours for the parameters / and 
Zc displayed in Fig.(j4]). It is clear that adding the BAO 
points enables us to increase slightly the precision in 
the determination of the free parameters. 

A final comment would be useful here. A direct com¬ 
parison between Fig. and Fig.® is not completely 
fair. In fact, the LOSS sample we have used, does not 
include the estimated systematic errors, whereas the 
JLA set does. This is the reason for the difference be¬ 
tween the precision in the parameter determination ob¬ 
served in the graphs. 
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f 


Fig. 3 We plot the contour at la and 2a for the param¬ 
eters / and 2c after marginalizat ion of M. Here we have 
used data only from the JLA set (iBetoule et aL|[20l4l . 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 


f 

Fig. 4 We plot the contour at la and 2a for the param¬ 
eters / and Zc after marginali zation of M. Here we have 
used SNIa data from the JLA (jBetoule et ahllioidl l and the 
BAO points. 


4 Discussion 


From a theoretical point of view, we know that the in¬ 
termediate evolution is more appropriated to describe 
the low redshift evolution, because this allows the ex¬ 
istence of a transition from the decelerated expansion 
phase to the current accelerated expansion phase. 

In previous sections, we showed that the intermedi¬ 
ate evolution successfully fit the data from SNIa and 
BAO. In fact, for the best fit using SNIa (JLA set) 
only, we get Xmin = 32.89 for the intermediate evolu¬ 
tion (with three free parameters, /, Zc and M), whereas 
within the ACDM (with also three free parameters (0^, 
Oa and h j) we get Xmin = 32.76, both certainly a good 
fit to the data (for the binned JLA set with 31 points, 
the estimated variance in is cr ~ a/ 2/31 ~ 0.25, 
whi ch means that both fi t are comparable to each other 
(see Andrae et al. ( 2011)11 11. 

Just to make the comparison more evident, we plo t 
the binned data from the JLA set ( Betoule et al.ll2014ll , 


with the best fit curve obtained from the theoretical 
model we studied - the intermediate evolution Eq. (HD 
- and the po wer-law model a(t) ~ with (3 = 1.55 
(quoted from ( Dolgov et al. 2014ll l in Fig.([5|). 



Fig- 5 We plot the 31 data points from the JLA set 
dBetoule et al.ll2014l l together with the best-fit value found 
in this work, using the intermediate e volution (fTIl (^contin - 
uous line), and the best fit found in (iDolgov et aLll20l4 l 
using a power-law evolution a{t) ~ with = 1.55 quoted 
from that paper (dashed line). Notice how the power-law 
fit moves away the data points from 2 > 0.2. 


It is clear that the intermediate evolution provides a 
better fit to the data, than the power-law with /? = 1.55. 

Actually, our study suggests that if we only use SNIa 
data, using all the data from the LOSS compiled set, 
and independently using all the data from the com¬ 
pressed JLA set, the power-law best fit does not corre¬ 
spond to /3 ~ 3/2 rather we obtain /3 = 1.78±0.09 with 
LOSS, whereas /3 = —0.35 ± 0.06 with JLA. 

In this paper we have studied the intermediate evolu¬ 
tion as the most appropriate for describing the low red- 
shift regime supported by observations from type la su¬ 
pernovae and BAO. We found th at the recent data (us - 
ing the LOS S compiled sainple (^Ganeshalingaml 120131) 
and JLA set ( Betoule et al.ll2014ll l suggest that the in¬ 
termediate evolution characterized by the scale factor 
(HD is as good a fit to this redshift regime as the stan¬ 
dard ACDM model. 

This result can be considered a step forward from 
( Dolgov et al. 2014[) . where the authors used a power- 
law type evolution a{t) ~ and studied the constraints 
the data imposed on the parameter (3. Although they 
successfully produce a good fit to the SNIa data, the 
model they consider is unable to describe a transition 
from a decelerated to an accelerated expansion phase. 
This transition is a key ingredient that any model must 
satisfy if it claims to describe the recent (low redshift 
2 < 1) evolution of the universe, which is essentially 
the same redshift range as the supernova data span. 

Although these results seems to rule out the power- 
law evolution to describe the low redshift regime (unless 















































6 


we are interested in fitting the very low regime 2 < 0.2), 
the original motivation, that of confronting modified 
gravity theories to data looking for signals in favor of a 
dynamic adjustment mechanism to the vacuum energy 
problem, remains intact. Actually, this becomes even 
more interesting with the connection with the Gauss- 
Bonnet coupling terms needed in the action to get the 
intermediate evolution. 
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